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Abstract 

Future electron-ion colliders (elC) will focus on the unitarity properties of deep inelastic 
scattering (DIS) in the limit of strong nuclear absorption. Strong nuclear shadowing and 
a large abundance of coherent diffraction are the most striking consequences of unitarity, 
and here we report quantitative predictions for these effects in the kinematical range of the 
planned elC. 
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1. Introduction. 



In deep-inelastic scattering off nuclei at the Bjorken variable 



1 



= 0.15A- 1 / 3 , 



(1) 



x < x A = 



m N R A 



where Ra is the radius of the target nucleus of mass number A and is the nucleus mass, 
unitarity driven nuclear shadowing (NS) becomes important [1, 2, 3]. It comes along with the 
large-rapidity-gap coherent diffractive DIS (LRG DIS) when the target nucleus remains in 
the ground state. For strongly absorptive nuclei the unitarity condition was shown to imply 
a paradoxically large, ~ 50%, fraction of LRG DIS [4] (for the experimental confirmation of 
LRG DIS off nuclei see [5], similar calculations were reported in the recent Ref. [6]). It is 
precisely the s-channel unitarity condition which controls the interplay of the virtual and real 
pQCD radiative corrections in the small-x evolution of different nuclear observables [7]. In 
this communication we report quantitative predictions for the small- a; evolution of NS and 
LRG DIS for the energy range of future elC [8] which will test the unitarity properties of 
hard scattering processes under strong nuclear absorption (saturation) . 

2. First iteration of the LL(l/x) evolution for nuclear cross sections. 

We base our work on the color dipole (CD) approach to the Leading Log(l/x) (LL(l/x)), 
or BFKL [9], evolution of DIS [10, 11]. For nuclear targets a complete resummation of LL(l/x) 
effects is as yet lacking, the Double Leading Log [12], large- N c Color Glass Condensate [2] 
and the fan diagram resummation [13] approximations were considered in the literature ( for 
the review see [14], the dominance of fan diagrams was questioned, though [15]). Our point 
is that at energies of the planned elC [8], the x-dependence of NS and LRG DIS is practically 
exhausted by the first CD LL(l/x) iteration which is calculable exactly without invoking the 
large-iV c approximation [2, 3]. Indeed, the first CD LL(l/x) iteration dominates in the x 
region 
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Figure 1: The predictions form CD LL(l/x) evolution vs. the NMC data on the shadowing ratio 
Rsn/c [26]. The solid lines represent the full shadowing from the qq and qqg Fock states, while the 
dashed lines show a shadowing at the qq level, see also a caption to Fig. 3. 

where A e // ~ 0.1 — 0.2 is the exponent of the local x-dependence of the proton structure 



function, F 2p oc x 



-A 



[16]. 



When viewed in the laboratory frame, NS derives from the coherent interaction of qq, qqg, ... 
states. To the required accuracy, the Fock state expansion of the physical photon [7*) reads 
|7*) = \J~Z~ g ^ qq\qq) + & qq - g \qqg), where fy qq - and § q q g are the light-cone wave functions (WF's) 
of the qq and qqg states, \J~Z~ g is the renormalization of the qq state by the virtual radiative 
corrections for the qqg state. For soft gluons the 3-parton WF takes the factorized form 



§ qq -g = ^qq{^qg — ^ qg} [10, 11]. The nuclear coherency condition reads [1] 



x/P <x A , 



(3) 



where x = Q 2 jlm^v is the Bjorken variable, v is the energy of the photon, M is the invariant 
mass of the multiparton Fock state and (3 = Q 2 /(M 2 + Q 2 ). 

In the CD BFKL-Regge phenomenology of DIS one usually formulates the boundary 
condition at x = x = 0.03 [17]. For extremely heavy nuclei xa <^ x and the LL(l/x) 
evolution starts at x = Xa with the boundary condition formulated in terms of the free- 
nucleon quantities LL(l/x)-evolved from xo down to x A . At small x satisfying the condition 
(3) the color dipoles {r n } in the multiparton state are conserved in the interaction process. 



3 



R, 



0.6 



I I I mini — I I I I I mi — I I I I llll 



C/D 




I I I I 



~l 1 1 — I — llll 



Ca/C 




4 -3 

10 10 



10" 



10" 



J I I I llll 



10" 



Figure 2: The same as in Fig. 1 for the NMC data on the shadowing ratio Rc/d [27] and Rca/c 
[28] as a function of two correlated variables x and Q 2 (for details of x — Q 2 correspondence see 
[27], [28]). 



At the boundary x = x A and for A ^> 1, the nuclear S-matrices equal [18, 19] 

S n (x A , b, {r} n ) = exp[-^a n (x A , {r} n )T(b)], 



(4) 



where a n (x A , {r}„) is the free-nucleon CD cross section for the n-parton state, and [20] 

ct 2 ,a^a, r) = 2 J d 2 b[l - S 2 (x A , b, r)]. (5) 



Here T(b) = dzriA(Vz 2 + b 2 ) , is the optical thickness of a nucleus at an impact parameter 
b, the nuclear matter density n^(r) is normalized as / d 3 rn(r) = A. 

A quite common extension of Eq.(5) to x x A with the BFKL-evolved a 2 (x,r) is com- 
pletely unwarranted. Instead one must evolve the nuclear S-matrix. Specifically, after the 
gluon variables have been properly integrated out, the effect of the extra gluon in the Fock 
state expansion for the incident photon boils down to the renormalization of the S-matrix 
and nuclear cross section for the qq CD [10, 11, 7],: 

dS2{X d p h) = Jd 2 pmp)-^(p + r)\ 2 [S 3 (x,p,r;b)-S 2 (x,r;b) 
a A (x,r) = a 2 ,A(xA,r) + a A \x, r), 
v$\x,r) = 2\og{^jd 2 bjd 2 p\^{p)-^p + T)\ 2 
S 2 {x A) r; b) - S 3 (x A , p, r; b) , 



x 



N 2 



N 2 -l 



a 2 (x, p) + a 2 (x,p + r) 
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N 2 -l 



a 2 (x,r) 



(6) 



where p is the qg dipole and 



y/C F a s p p . p , , 



is the radial WF of the qg state with the Debye screening of infrared gluons [10, 11]. The 
virtual photoabsorption cross section is an expectation value <7a(x, Q 2 ) = {qq\a A (x,r)\qq) . 
The shadowing ratio Ra(x,Q 2 ), decomposed into the qq and qqg contributions, equals 

D/ ^ cr A (x,Q 2 ) Aa 2tN (x A ,Q 2 ) - a 2A (x A ,Q 2 ) 

Ka{X, Q ) = — — = 1 — 



Aa N (x,Q 2 ) ' Aa N (x,Q 2 ) 

, (xa\ Aa^\x A ,Q 2 )-a A l \x A ,Q 2 ) 
" l0g \~) Aa N {x,Q 2 ) ' (8) 

where a N (x,Q 2 ) = <J2,n(xa,Q 2 ) + ^og(x A /x)a ( ^ > (xa, Q 2 ) is the LL(l/x)-evolved free-nucleon 

cross section. Our interest is in the well evolved shadowing at i < x^, the onset of NS at 

x > xa must be treated within the light-cone Green function technique [12, 21], here we show 

only the gross features of the large-x suppression of NS following the prescriptions from Refs. 

[20, 22] and the /9-dependence of LRG DIS as predicted in [23] and confirmed in the HERA 

experiments [24]. We use the free-nucleon CD cross section tested against the experimental 

data from HERA, it is described in the Appendix. The nuclear density parameters are taken 

from the compilation [25]. In Fig. 1 we compare the predictions from the first LL(l/x) 

iteration with the NMC data [26] on the ratio of R Sn /c(x, Q 2 ) = R Sn (x, Q 2 )/R c (x, Q 2 ). The 

results for Rc/d [27] and Rca/c [28] are shown in Fig. 2. The agreement with the experimental 

data is good, but in a limited region of x, Q 2 the contribution from the shadowing of the qqg 

Fock states is still small. A much broader range of x and Q 2 can be covered at elC [8], and 

in Fig. 3 we show our predictions for Ra(x, Q 2 ) at elC. 

3. LL(l/ar) evolution for LRG DIS 

The LL(l/x) evolution of fully inclusive forward LRG DIS off a free nucleon, starting from 
xjp = x , has the Fock state expansion [23, 11] 

da®(x,x ) 



(qq\cr 2 (x )\qq) + (qqg\cr 3 (x ) - (T 2 (x )\qqg) + ... , (9) 



dt 

where t is the (p,p') momentum transfer squared. Upon the integration over the gluon mo- 
mentum the qqg term gives rise to the LL(l/x) evolution of LRG DIS. Its decomposition 




-4 -3 -2 -4 -3 -2 -4 -3 -2 -1 

10 10 10 10 10 10 10 10 10 10 

X 

Figure 3: The upper panels: predictions from CD LL(l/x) evolution for NS ratio Ra for C, Ca and 
Pb nuclei as a function of x at Q 2 = 1, 5, 20 GeV 2 , the bottom to top curves, respectively. The solid 
lines show the combined NS (8) from the qq and qqg Fock states of the photon, while the dashed 
lines represent NS at the qq level, i.e., the term i?2,A = 1 — [A(T2,n(xa, Q 2 ) — <?2,a{xa, Q 2 )]/ 'A(jn{x) 
in the expansion (8). The lower panels: the same as above for the fraction of DIS which is LRG 
coherent diffraction, R co h- Dashed lines represent the contribution from the low-mass qq excitations, 
(?a M , the sum of the low-mass and 3IP high- mass terms is shown by solid lines. 

[29] 

of (x ) - 02(2:0) = [03(2:0) - 02(a; o )] 2 + 2a 2 (x ) (03(^0) - 02(^0)) (10) 

splits the LL(l/x) evolution into the real production of high-mass (triple-pomeron, 3D 3 ) 
states governed by [os(xo) — o- 2 (xo)] 2 and the virtual pQCD radiative correction to the Born 
excitation of low-mass, M 2 ~ Q 2 , states described by the first term of eq.(10). 

For heavy nuclei we consider coherent diffractive DIS per unit area in the impact parameter 
space, da® (x, xq) / d 2 b with the boundary xq = xa- At fixed x <^ xa, the small-mass qq diffrac- 
tion enters at xip ~ x, while the small-x evolution of the high-mass diffraction starts with the 
rapidity gap variable xjp = xa- To the LL(l/x) approximation, integrating out the gluon vari- 
ables, one can cast it in the color dipole form [29] da^(x, xa, b)/d 2 b = (qq\T,^(x, xa, b, r)\qq), 
where, with certain reservations [30], one can dub S^(x,XA,b,r) a CD cross section for the 
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Figure 4: The predicted nuclear mass-number dependence of the ratio, R(3JP /LM), of the high- 
mass to low-mass diffractive DIS at several values of Q 2 . 

pomeron target. Extending the decomposition (10) to nuclear targets, for the high-mass 
excitations, i.e., 3IP at x <C x A , one finds [29] 



(qqg\ [S 2 (x A ,b,r) - S 3 (x A , b, p, r)} z \qqg) = (qq\E A w (x,x A ,b,r)\qq) 



Ef{x,x A ,b,r) = log(^) /rf 2 p|^(p)-^(p + r)| 2 [S 2 (xA,b,r)-S 3 (x A ,b,p,r)] 2 (li; 
Analogously, for the low-mass (LM) excitations da L A M ' jd 2 b = (qq\^ A M \qq) and 



ET(x,x A ,b,r) = [l-S 2 (x A ,b,r)]^ + 21og(^) [1 - S 2 (z A , b,r)] 



x A 



X 



J d 2 p\i,{p) - V(p + r) | 2 [S 2 (x A , b, r) - S 3 (x A , b, p, r)] . (12) 



A fraction of DIS which is coherent diffraction, R co h(x,Q 2 ) = a A (x,x A ,Q 2 )/cr A (x,x A ,Q 2 ), 
is shown in Fig.3. The large- X suppression of R co h(x,Q 2 ) by the nuclear form factor is eval- 
uated with the mass spectrum [23, 24]. For the realistic dipole cross section even the lead 
nucleus is still a grey one: the predicted R co h is substantially smaller than the black disc 
result R coh = 0.5 [4]. 

The total cross section of the low-mass diffraction a A M (x, Q 2 ), i.e., DIS off the large-/3 
valence qq state of the pomeron, is saturated by the contribution from M 2 ~ Q 2 . The strength 
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of 3IP diffraction, or DIS off the small-/? sea in the pomeron, is conveniently measured by 

GwM = |/ <Pb(q$HF(x,x A ,b,T)\q$ = (M 2 + Q 2 ) da ^ X 2 F \ (13) 

In Fig. 4 we show the ratio R(3iP/LM) = G 3 jp(xjp)/a^ M (x), where we take x = 10" 3 and 
xjp — x/P — 10~ 2 , i.e., (3 = 0.1, where the diffractive excitation of the qqg states takes over 
the excitation of the qq states ([23], for the experimental data see [24]). At such a value of 
xjp the nuclear form factor effects can be neglected. The rise of R(3TP/LM) with Q 2 is a 
standard pQCD scaling violation growth of the sea parton density. Subtleties of the scaling 
violations at small Q 2 close to the nuclear saturation scale Q\ will be discussed elsewhere. 

The suppression of R{?>W / LM) for heavy nuclei, evident in Fig. 4, is a unitarity effect. In- 
deed, the nuclear absorption in the integrand of triple-pomeron cross section, H\ w (x, xa, b, r), 
is stronger than that of S^ M (x, xa, b, r), and it is steadily enhanced for large A. 

5. Summary. 

Based on the color dipole approach which correctly reproduces the experimental data on 
the proton structure function measured at HERA and the NMC data on nuclear shadowing, 
we reported the quantitative predictions for the LL(l/x) evolution of the nuclear shadowing 
and coherent diffraction dissociation off nuclei in the kinematical range to be covered at future 
electron-ion colliders. We explicitly separated the contributions to both observables from the 
low-mass and the triple-pomeron high-mass diffractive states, including the virtual radiative 
correction to the low-mass diffraction. Regarding the saturation properties, even the lead 
nucleus is still a grey one: the fraction of DIS which is diffractive is substantially below the 
black disc result R co h = 0.5. The gross features of the predicted x and Q 2 dependence of 
nuclear shadowing do not change much from the carbon to lead target; we also predict a 
suppression of the high-mass vs. low-mass coherent diffractive DIS on heavy nuclei. 

Acknowledgments. This work has been partly supported by the DFG grant 436 RUS 
17/82/06. VRZ acknowledges also partial support from the RFBR grant 06-02-16905-a. 

Appendix. Boundary condition for the LL(l/x) evolution. 
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Figure 5: The proton structure function F 2p (x,Q 2 ) predicted by the BFKL-evolved color dipolc 
cross section (14) vs. the experimental data from ZEUS (open circles, [31]), HI (triangles, [32]) and 
E665 (squares, [33]) Collaborations. 



The realistic input a(xo, r) must be denned for all scales r, from small to non-perturbative 
large. Motivated by our successful BFKL-Regge color-dipole phenomenology of DIS [17, 29], 
at xq = 0.03 and in the limited range of r > 10~ 3 fm of the practical interest, we take 

where = (r/r ) 2+7 % r = 1 fm, a 1 = 22.5 mb, a 2 = 20 mb, a 3 = 8 mb, 71 = 0.15,7 2 = 0.8, 
73 = 2. Its LL(l/x) evolution is described by the Color Dipole BFKL equation [10, 11] with 
the infrared freezing of the one-loop, 3-flavor QCD coupling, 

®s(q 2 ) = - r,. 47r o,,.o T ' ( 15 ) 



91og[(g 2 + g|)/A| CD 

where qf ~ 0.7 GeV and Aqcd — 0.3 GeV. A new observation is that at the expense of 
somewhat smaller Debye screening radius, R c = 0.2 fm vs. R c = 0.275 fm in [10, 11], the 
LL(l/x) evolution of the boundary condition (14) gives a good description of the experimental 
data from HERA (see Fig 5) without splitting the boundary condition into the BFKL-evolving 
perturbative, and non-evolving, non-perturbative, components. 
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